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Abstract 

For non-critical string theory the partition function reduces to an integral over moduli 
Ph space after integrating over matter fields. The moduli integrand is known analytically 

for genus one surfaces. The formalism of dynamical triangulations provides us with 

a regularization of non-critical string theory and we show that even for very small 
t— H triangulations it reproduces very well the continuum integrand when the central charge 

c of the matter fields is large negative, thus providing a striking example of how the 
(V) quantum fluctuations of geometry disappear when c — > — oo. 
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1 Introduction 



Non-critical string theory (or equivalently 2d Euclidean quantum gravity coupled to 
conformal field theories) has allowed us to study non-perturbative aspects of string 
theory and quantum gravity. Certain aspects of the theory can be solved by continuum 
methods and certain aspects have been solved using combinatorial methods (matrix 
models and dynamical triangulations (DT)). Originally dynamical triangulations were 
considered as a lattice regularization of path integral over two dimensional world sheet 
geometries [H [2J, [3] , but the versatility of the matrix models allowed for many other 
interesting interpretations, which became increasingly detached from the original idea 
of DT as a lattice regularization of a 2d quantum world. 

A typical dynamical triangulation represents a (piecewise linear) geometry which 
appears in the path integral over 2d geometries. As such it does not represent anything 
like a physical measurable geometry, in the same way as a random path in the path 
integral representation of the particle is not physical and cannot be measured. However, 
it has interesting fractal properties, precise as the random path for the particle with 
probability 1 has Hausdorff dimension 2. The "random geometries" in the path integral 
(represented in DT by dynamical triangulations) have Hausdorff dimension 



Contrary to the particle case, this is an intrinsic fractal dimension, not a fractal dimen- 
sion in target space. Also, it depends on the central charge of conformal field theory 
coupled to 2d gravity. This result was first obtained for c = in [3] using the DT 
formalism, and later, again using the DT formalism, from a different perspective in 
[5j E]. The general formula was derived in [7j using Liouville theory. 

These results indicate that for a given central charge c there exists a measure defined 
on the set of continuous geometries. However, no precise mathematical definition exists 
yet, contrary to the situation for continuous paths (the Wiener measure). In critical 
string theory this problem is circumvented by using the conformal invariance of the 
theory. In this way most of the integration over geometries can be factored out, and 
we are left only with a finite dimensional integration over the moduli parameters (and 
a non-trivial integrand for higher genus surfaces). In non-critical string theory the 
situation is more complicated due to the zero mode of the Liouville field, except in the 
genus one case. In this case one obtains just the standard result from critical string 
theory [8]. Thus, if we consider the situation where we fix the world sheet area to be 
A we obtain 

Z^ h=1 \A) ~ A^ 1 [ ^F(r) c -\ (2) 
where r = T\ + ir 2 is the moduli parameter, Ai denotes the fundamental domain of r 
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Figure 1: The fundamental domain Ai in the upper-half complex plane. Its colloring 
illustrates the choice of bins used in Fig. [3j The maximum of F{r) l ~ c at r max = e l7r//3 
corresponds to the equilateral torus, i.e. the quotient of the plane by the hexagonal 
lattice. 
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The integral in ^ diverges for c > 1, which tells us that there is a c = 1 barrier 
beyond which the formulas make little sense. When c decreases geometries with large 
T2 become more and more suppressed. In fact the maximum of F(r) 1_c , located at 
e W3 j n the fundamental domain, will dominate the integral for c — > — oo (see 
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Fig.0. 

By assigning moduli parameters to triangulations of the torus one can study how 
the formalism of DT reproduces the integrand in (|2]). The possibility to assign a 
moduli parameter to a toroidal triangulation was first realized in [9], where the first 
study the DT-realization of eq. ^ was performed. Using the formalism of simplicial 
cohomology this assignment can be done in a natural and relatively simple way [TO] , 
and one can even explicitly construct the harmonic map from the geometry defined 
by the triangulation to the parallelogram defined by the moduli parameter. Stated 
briefly, given two independent cycles on a torus, 71 and 72, one can find dual harmonic 
one-forms ati and a 2 such that 

' « J = ^ (5) 
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1 A convenient choice for the fundamental domain A4 is (see Fig. [TJ): 



r e M if 



r 2 > 0, - i < n < and |r| > 1 
r 2 > 0, < ti < i and It > 1 



(3) 
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Then a modular parameter of the torus is given by 

T ~~ (a 2 |a 2 ) V (a 2 |a 2 ) V(" 2 |« 2 )/ ' 

where (•, •) refers to the standard inner product on one- forms. This r is not necessarily 
in the fundamental domain A4, but is uniquely associated with a point there by a 
modular transformation. The concepts of cycles, one-forms and harmonic one-forms 
have a natural formulation on complexes, the de Rham cohomology being replaced 
by simplicial cohomology. In the case of DT, where the triangles are equilateral, all 
formulas become very simple. On toroidal triangulations we can first find cycles (non- 
contractible link- loops), then the dual one-forms (harmonic in the simplicial context), 
the corresponding r and finally the modular transformation which maps it to the 
fundamental domain. We refer to [10] for details. 

In [TQl |TT] this formalism was used to study 2d Euclidean quantum gravity for 
the matter central charges c = and c = — 2. These values were chosen since they 
allowed for large computer generated triangulations with the correct weights. Large 
triangulations were necessary in order to properly include the singular geometries with 
large t%. In this article we extend the study to large negative c. For such c we expect 
geometries with large T2 to be suppressed and therefore smaller triangulations should 
suffice. As we will explain below, relatively small triangulations can be computer 
generated with the correct weight for any c. 



2 Large negative valued central charge 

Consider the partition function for the bosonic string in d dimensions: 

Z™ = J V[g] J V g X,e~ s ^\ (7) 

where the integration is over worldsheet geometries [g] of genus h and Gaussian matter 
fields X^. The d Gaussian fields can be integrated out and we obtain 

|p 9 X,e-™~(det(-A;))^ /2 . (8) 

Here denotes the Laplace-Beltrami operator on the background geometry defined 
by a the metric g a b and the prime signifies that the zero mode has been removed when 
calculating the determinant of A g . If we consider d as a formal parameter we can 
continue this expression to negative d. For positive integer d we can identify d with 
the central charge c of a conformal field theory coupled to the 2d world sheet geometry 
and formally we can do the same for any real d. The DT-formalism tells us that we 
should represent the regularized partition function (J7|) as 

Z^\N) = J2?^- (det(-A' TN ))' d/2 , (9) 



4 



c = -5 c = -10 c = -20 c = -40 c = -80 




Figure 2: Typical triangulations appearing in the partition function with N = 160 
triangles and increasingly negative central charge c. Thick (red) curves correspond to 
the shortest non-contractible loops studied in section [3] 

where the summation is over all triangulations of the torus with iV triangles. CV is 
the symmetry factor of the triangulation T, i.e. the order of the automorphism group 
of the graph T. Finally Ay denotes the (discretely defined) Laplacian on the DT- 
surface, which we take to be the usual graph Laplacian of the 3 graph dual to the 
triangulation. The continuum area A in formula rt2J) is related to N by A = N a 2 \/3/2, 
where a denotes the length of a link in the triangulation constructed from equilateral 
triangles. The continuum limit of DT is obtained by taking a — > while keeping A 
fixed. Thus we only expect the DT result to be exact (i.e. to agree with the continuum 
expression) for N — > oo. 

As long as N is not much larger than a few hundred triangles standard Monte 
Carlo simulations can be used to generate an ensemble of triangulations {Tj(iV)} with 
the correct weight dictated by the partition function ([9]). We refer to the Appendix 
for details on the updating algorithm used in the simulations. For each triangulation 
Ti(N) we can now extract the moduli parameter by the technique mentioned above 
and described in detail in [TD]. As mentioned before one can obtain an explicit har- 
monic map from the triangulation to the parallelogram in the Euclidean plane defined 
by the moduli parameters. Fig. [2] shows examples of such maps for typical triangu- 
lations appearing in the partition functions for various valus of the central charge. 
These examples indicate that typical triangulations become more and more regular, 
i.e. closer to the hexagonal lattice in Fig. [TJ when c becomes more negative. From the 
measurements we construct a probability distribution in the fundamental domain A4 
which we can compare with the continuum integrand F(r) c_1 /r| in The errors 
one encounters are twofold. There will be statistical errors which diminish with the 
size of the computer generated ensemble {T/v} in the standard way, and there will be 
systematic errors associated with the use a finite number N of triangles. 

Our results are shown in Fig. [3} The figures are constructed in the following way. 
The fundamental domain Ai is partitioned into bins according to the value of T2 with 
bin size AT2 = 0.025, as illustrated in Fig. [I] (however notice that there we have 
displayed bins with larger bin size for visualization purposes). We have chosen to 
disregard the value of T\ since the theoretical integrand F(r) c_1 /r| varies little with 
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Figure 3: The r 2 distributions rf p(r 2 ; N) compared to the theoretical distributions. 
The data points and the theoretical curves are produced as described in the main text. 



T\. We denote by p(r 2 )Ar 2 the probability of observing a r(Tj(iV")) associated with a 
triangulation Tj(iV) sitting in the bin r 2 — ^Ar 2 < T 2 (Ti(N)) < r 2 + |At 2 . We compare 
p(r 2 ) to the theoretical distribution Ptheory^) obtained by integrating F{r) c ~ l /rf over 
the corresponding bin. The measured distributions for various values of the central 
charge are depicted in Fig. [3] by points with error-bars (which are hardly visible for 
c < —10) and the solid lines correspond to the theoretical distributions. 

From Fig. [3] it is seen that r|p(r2, N) is more and more concentrated at the smallest 
possible r 2 values when c becomes increasingly negative, as expected from the theo- 
retical distribution. Also, in order to obtain a reasonable accuracy we have to take 
iV = 300 for c = —5, while for c = —40 one obtains good results already with iV = 60. 
The need for large N the closer c is to barrier c = 1 reflects that the distribution 
p(r 2 ) becomes more flat when c increases. Thus both (r 2 ) and (rf) — (r 2 ) 2 increases. 
However, the large r 2 region is not well reproduced by triangulations with small N, the 
reason being that loops take integer length and, loosely speaking, r 2 is related to the 
ratio between the two independent shortest non-contractible loops. When we get closer 
to c = 1, the effective area of the fundamental domain M. that has to be covered by 
the simulations becomes larger. Because of this and the need for larger triangulations, 
more and more computer-time is needed. This trend is reflected well in the figures and 
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the message is: the semi-classical region c — > — oo is easy! 

As further evidence of the precision with which the simplicial cohomology is able to 
reproduce the continuum results in the range of central charges that we have considered 
here, we have fit the data in Fig. |3]to the theoretical distribution with the central charge 
c in formula (]2]) as a free parameter. The result is listed in Table [I] 



c 


-5 


-10 


-20 


-40 


Cfit 


-5.01 ±0.1 


-9.94 ±0.04 


-19.9 ±0.1 


-39.88 ±0.05 



Table 1: The fit cju to the integrand in pi) from data. 



3 The Hausdorff dimension 

One of the most remarkable formulas in 2d Euclidean quantum gravity is Why do 
we need to verify it? There have been a number of competing formulas [13]. Although 
it is mostly believed today that these formulas describe critical properties of the matter 
theories coupled to 2d gravity rather than the fractal structure of a typical geometry 
which appears in the path integral [H], the question is not completely settled. In 
addition, the derivation of ([!]) is based on a certain number of assumptions related to 
diffusion on fluctuating geometries, assumptions which are not necessarily valid. 

In pH CLE] formula Q was verified with great numerical accuracy for c = and 
c = — 2. Here we check the formula for more negative c. As noted in [T0J E] the torus 
offers a nice opportunity for a numerical check of the formula which is independent of 
the original numerical check [T2] performed by simply measuring the area A(r) enclosed 
in disks of geodesic radius r: 

(A(r))~ r dh , r^A 1/d \ (10) 

where A is the fixed area of the 2d geometry, as in formula (J2]). On the torus one can 
use that in a given geometry a shortest non-contractible loop must be a geodesic curve. 
On average a non-contractible loop will have a length (L)a ~ A 1 / 2 . However, if the 



geodesic distances scale anomalously, as is the case if dh ^ 2 in ( 10), we expect instead 



(L)a ~ A 1/dh . (11) 

In DT we can perform such measurements: we define curves as link-paths, and a 
geodesic curve between two vertices as a shortest link-path. Let now {T;(iV)} denote 
the DT ensemble of triangulations with N triangles generated with the correct weight. 
Fixed N corresponds to fixed continuum area A in this context. For each triangulation 
we determine a shortest non-contractible loopj^] and in this way we can calculated the 

2 The practical algorithm used to find the shortest non-contractible loop uses in a nice way simplicial 
cohomology. We refer to [10] for details. See also Fig. [2] for examples of shortest non-contractible 
loops. 
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Figure 4: Left figure: a log-log plot of ((L — a)) N versus N as well as the theoretical 
curve bN l l dh for c = —20, a = —0.10, b = 0.61 determined by a best fit. The range of 
N is 30-400. Right figure: the difference ((L - a)) N - bN l ' dh . 



average (L)n. We do this for as large range of iV as possible and determine dh from a 
fit to (11). We have chosen c = —20 and N in the range from 30 to 400. A best three 
parameter fitj^] 

((L - a)) N = b N l/d d = 2.76 ±0.07, (12) 

which should be compared to the theoretical value dh ~ 2.66. In view of the small 
values of iV used we consider this result as quite satisfactory and a confirmation of 
formula Q. The left figure in Fig. [4] shows the measured ((L — a))jv together with the 
theoretical curve bN l l dh (where a, b is determined by a best fit). The right figure in 
Fig. [4] displays the error in the fit. It is clear that N is too small to produce a perfect 
fit. If one leaves out the three smallest values of iV one obtains a \ 2 value of order 2, 
which is good. 



4 Discussion 

Matrix models have provided convincing evidence that the formalism of dynamical tri- 
angulations provides a regularization of non-critical string theory or equivalently 2d 
Euclidean quantum gravity coupled to matter with central charge c < 1. However, 
quantum geometric aspects of 2d Euclidean quantum gravity are not easily analyzed 
using matrix models. The reason is that when using matrix models one integrates over 
all geometries. The matrix model thus provides us with the full partition function of 
matter and gravity and to obtain more detailed information one has to insert certain 
"punctures" and calculate their expectation values. These operators also mix matter 
and geometry and to disentangle geometry and matter is not easy and becomes in- 
creasing complicated in the limit c — > — oo, a limit where one naively would expect 



The parameter a in (12) is not present in (11), but is quite natural. It reflects the fact that L 



in DT is an integer, contrary to the L in the continuum formula. It can be viewed as a finite size 
correction, using the terminology of finite size scaling in the theory of critical phenomena. 
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Figure 5: Flip move. 



the quantum aspects of geometry to decouple. Formula ([I]) can be read as an explicit, 
quantitative statement about this decoupling: d h — > 2 for c — > — oo. However, it is not 
derived using matrix models but quantum Liouville theory. 

One would expect that the original interpretation of DT, where the triangulations 
are viewed as a lattice regularization of the path integral over geometries, should allow 
us to observe the increased semi-classical nature of geometry in detail. The purpose 
of this article was to provide evidence for this. For the torus we have shown that the 
formalism of DT indeed reproduces the integrand in ([2]), not only the partition function. 
Further, the more negative the value of c, the better is the continuum expression 
reproduced even for small triangulations. While conformal invariance of the continuum 
theory blurs somewhat the concept of a regular geometry, the possibility to use the 
moduli parameter of the torus has allowed us to address the problem in a conformally 
invariant way, and in the c — > — oo limit we clearly see the dominance of the "regular" 
triangulations, i.e. triangulations of the torus which have small values. 

Finally, using the concept of a shortest non-contractible loop on the torus, we have 
been able to test and (to some degree) verify formula Q. 

Acknowledgment 

The authors acknowledge support from the ERC-Advance grant 291092, "Exploring 
the Quantum Universe" (EQU). JA acknowledges support of FNU, the Free Danish 
Research Council, from the grant "quantum gravity and the role of black holes" . Part 
of this research was carried out at Utrecht University, The Netherlands. 

Appendix: Monte Carlo simulations 

In order to generate random triangulations according to the distribution ^ we use 
standard Monte Carlo methods. We start with some initial triangulation with the 
desired number N of triangles. Then we apply random updates to the triangulation via 
so-called flip moves, as shown in Fig. [5j To obtain the correct Boltzmann distribution 
one has to implement an acceptance probability that is related to the ratio of Boltzmann 
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factors before and after the flip move, i.e. 



V 



det'(-A - 6 A) 
det'(-A) 



-d/2 



[det(/ + A^A 



-d/2 



(13) 



Here 5 A is the change under the flip move of the Laplacian on the dual graph of the 
triangulation, and A -1 is the pseudo-inverse of A. The matrix 5A has only eight 
non-vanishing entries and, adopting the labeling in Fig. |5j is given by 



5A 



( 



-1 
1 

V : 






1 

-1 



-1 
1 






1 

-1 






Therefore (13) reduces to a determinant of a four by four matrix 

-d/2 

1 



V 



1>1 

'03 



V>1 
"^3 



where the vectors <fi and ip are the solutions to the linear equations 



A0 



/ 1\ 

-1 



V 7 



Ail) 



( 0\ 



1 

-1 

V '■} 



(14) 



(15) 



(16) 



Solving these equations using numerical linear algebra techniques is the most time- 
consuming step of the update move. It limits the simulations to triangulations of up 
to several hundred triangles. 
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